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The lattice implementation of monopole creation operator proves to have problems related to bulk 
transitions that can possibly affect the interpretation of its mean value as an order parameter for 
monopole condensation. Preliminary evidence is presented that these unexpected behaviours are 
in fact only due to lattice artefacts and do not spoil the physical interpretation of the monopole 
operator. 
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1. Introduction 

The idea that color confinement can be produced by dual superconductivity of the vacuum was 
introduced in Refs. [1,2] and since then it attracted a big interest. However, despite its intuitive 
nature, this proposal is not easy to demonstrate, the main difficulty being that it is far from trivial 
to identify the magnetic degrees of freedom. 

In Ref. [3] the idea was advocated that any operator in the adjoint representation of the gauge 
group can be used as an effective Higgs field to identify the magnetic U(l) group (abelian pro- 
jection), physics being independent of that choice. This revealed to be true as far as monopole 
condensation is concerned while for monopole detection a strong dependence on the abelian pro- 
jection adopted is observed (see the discussion in [4]). 

The correct procedure to detect monopole condensation, that is Higgs breaking of the U(l) 
magnetic symmetry, is to evaluate the vacuum expectation value (vev) of a magnetically charged 
operator fx: this vev has to be zero in the normal phase and can be different from zero in the 
dual superconducting phase. In order to pursue this strategy a magnetically charged operator jx 
was introduced and developed in Refs. [5, 6, 7, 8, 9]. Similar constructions were developed in 
Refs. [10, 11]. 

The operator Ll, however, proved to a more rigorous analysis to have a bad infrared behaviour, 
i.e. a bad thermodynamic al limit, as shown in Refs. [12, 13]. 

In the following sections we will sketch the construction of the Ll operator and present a 
possible way to overcome these problems. 



2. Construction of ji 

The simplest way to construct a magnetically charged operator is to define it as the operator 
that adds a monopole to the state to which it is applied ([14]). This is easily done in the Schrodinger 
representation: if we denote by 11; (x) the canonical momenta conjugate to the transverse compo- 
nents of the gauge fields (x), the operator Ll defined by 

Li{y) =exp (ij d 3 xZ; ; (x,y)lI,(x)) (2.1) 

is the translation operator of the gauge field by b(x,y): 

M(y)|A(x)) = |A(x)+b(x,y)) (2.2) 

This is equivalent to say that the following commutation relations hold 

[A i (x),Li(y)]=b i (x,y)Li(y) (2.3) 
[n ; -(x),At(y)] = (2.4) 

Because of the linearity with respect to of the 't Hooft tensor, it follows from Eq. (2.3) that 

[Q,Ll(y)]=mLl(y) (2.5) 

where Q is the magnetic charge operator and m is the charge of the field b(x,y). The operator Ll is 
thus charged if b(x, y) is the field in x of a monopole located in y. 
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On the lattice the canonical momenta correspond (with the Wilson action) to the mixed spatial- 
temporal plaquettes and the operator in Eq. (2.1) can be rewritten (up to 0(a 2 ) lattice artefacts) as 

H = exp(-/3AS) A5 = £Tr{n ;o (n,0-n; o (n,0} (2.6) 

n 

where 

n,- (M) =C/ ; (n,?)C/o(n+f,/)t//(n ) ? + l) t f/ (n,O t (2.7) 

W iQ {n,t) =U i (n,t)U (n+i,t)M i (n + i,t)U i {r,t + iyu {n,t^ (2.8) 

M,(n,0 = ( eXp(/ ^ (n ' y)4>)lf? = ° (2.9) 
A ' [0 if t t^O 

and 4> is the generator of the gauge group which identifies the magnetic U(\). The operator of 
Eq. (2.9) corresponds to the choice of the so-called random abelian projection, which simplifies 
the numerical calculations ([9]). The fact that the operator in Eq. (2.6) adds a monopole to a given 
configuration can be shown directly by an explicit gauge transformation (see Refs. [7, 8]). If we 
denote by 5 the Wilson action, the expectation value of p is given by 

<">= jw^ (2 - 10) 

For computational reasons it is convenient to use, instead of (p), its logarithmic derivative 

p = -^log(n) = (S) s -(S + AS) s+AS (2.11) 

Indeed, since (p) = 1 at /3 = 0, the p vev can be reconstructed as (p) = exp y Jq p{x)Ax^ . 

The quantity p, which is the one measured in lattice simulations, should have a finite value in 
the thermodynamic al limit for temperatures smaller than the deconfinement temperature (T < T c ) 



SU(2) group 

Wu-Yang monopole of charge 4 

~< 1 1 1 1 1 1 r 




2.26 2.28 2.3 2.32 2.34 2.36 



P 



Figure 1: Value of the p observable defined in Eq. (2.1 1). 
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Figure 2: Bulk transition in G2 gauge theory 

and should develop a negative peak at T c , scaling with the appropriate critical indices (see e.g. 
Ref. [8]). For T > T c , if the dual superconductor picture is correct, (ju) = in the thermodynarnical 
limit and p — > — 00. 

The first simulations were done at A 7 , = 4 and different spatial extensions N s and seemed to 
agree with these expectations ([8, 9]). A more careful analysis performed with different values of 
AT, [12] showed that 

1 . the observed peak did not move by varying N t as expected from the known dependence of T 
on/3 

2. for T < T c the values of p, instead of tending to a finite limit as N s — > °°, keep growing 
negative with the volume, so that (ju) = also for T <T C 

3. a physical peak appears at the right position, superposed to the background of 1. and 2. 

The direct meaning of pL as an order parameter was thus spoiled. 

For abelian theories the operator in Eq. (2.6) clearly shows the correct behaviour for T < T c . 
Moreover it was shown in Ref. [7] to be equivalent to the one constructed in Ref. [15] by explicitly 
using the duality transformation. For non-abelian theories this is less clear: it is by now known 
([12, 13]) that the quantity p defined in Eq. (2.11) seems to diverge also for T < T c . 

This behaviour is shown in Fig. 1 for case of the group SU (2) on lattices of temporal extent 
N t = 4. The deconfmement transition is located, in the thermodynarnical limit, at fi c = 2.2986(6) 
and, as expected, p develops a negative peak in the neighbourhood of this j8 value (finite size 
corrections are obviously to be expected). However it looks like there is a strong volume dependent 
background also for T <T C , where p should have a well defined thermodynarnical limit. 

A plausible working hypothesis was that the source of this problem were the 0(a 2 ) lattice 
artefacts which are present in the definition Eq. (2.6). If this were true it seems more appropriate to 
use the Wu-Yang form of the abelian monopole for the field b(x,y) instead of the Dirac form used 
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Figure 3: Value of the p su b observable defined in Eq. (3.1). 



in the past. This is because in the Dirac expression 

(1+COS0' 



(2.12) 



rsinf? 

the field is large near the z axis also far away from the monopole, while in the Wu-Yang formulation 

:i-cos0) 

(2.13) 
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b s - 



rsinf? 

(1+COSfj) 

rsinf3 



0< <7i/2 + e 
n/2-e < 6 < n 



the only singularity of the field is at the monopole location. Indeed by using the field Eq. (2.13) 
the background is reduced typically by about a factor of 3. The general behaviour is nevertheless 
divergent in the thermodynamical limit, see Fig. 1. 



3. A possible solution 

If the divergence of p below j8 c is driven by lattice artefacts, it comes from short distances and 
therefore it has to be independent of the temperature, that is of the temporal extent of the lattice. 
This proves indeed to be the case, both with the Wu-Yang and with the Dirac monopole. 

As an example we show in Fig. 2 the values of p calculated for the G2 group on the lattices 
6 x 20 3 and 20 x 20 3 . The (first order) deconfinement transition of the model for N t = 6 is located 
between j3 = 1.39 and j3 = 1.4 (see Ref. [16, 17]), while the big bump near j8 = 1.355 is just the 
non-physical bulk transition (in fact it is just an analytic crossover). As expected for a purely UV 
effects the value of p does not depend on the temporal extent of the lattice. 

This suggests a possible method to obtain a reliable definition of p : subtract from the p values 
calculated on asymmetrical lattices (finite temperature) the values calculated at zero temperature: 

Psub = P - P(r=o) (3-1) 
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Figure 4: Rescaled values of p su t,. 

In terms of the pL operator this amounts to use ( i u)/(/i)r = o instead of just (ju). 

The results obtained by using the subtracted form Eq. (3.1) are shown in Fig. 3 and now the 
behaviour for j8 < j3 c is the correct one. Near the deconfmement transition it can be shown (see e.g. 
Ref. [8]) that the minimum of the dip in the p values scales as v , where Af s is the lattice spatial 
extent, while its width varies according to the scaling law N s '^ V . 

These scaling laws are shown to be satisfied in Fig. 4. The curves are drown by use of the 
known value of the deconfmement coupling, jS f = 2.2986(6), and v = 0.6301(4), which is the 
value of the v critical index for the 3d Ising universality class. For j8 < j3 c the scaling is very good 
near the transition, while for j3 > j3 f data are consistent with the expectation that p su b — > —°° as 
V — > °°. A careful analysis is on the way. 

4. Conclusions 

We have presented indications that the unexpected behaviour of the order parameter (ju) for 
monopole condensation in non-abelian gauge theories has to be ascribed to lattice artefacts. In 
order to disentangle these artefacts from the physical signal one possibility is to renormalize the 
expectation value of the jU operator by its zero temperature value. This procedure is shown to give 
reliable results, which nicely agree with the predictions based on universality arguments. Further 
studies on this line are on the way. 
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